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Abstract. We provide an overview of recent work exploring the quark-mass dependence of hadronic ob- 
servables and the associated role of chiral nonanalytic behavior due to the meson-cloud of hadrons. In 
particular, we address an issue of great current interest, namely the degree of model independence of 
results obtained through a controlled extrapolation of lattice QCD simulation results. Physical insights 
gained from this research are highlighted. We emphasize how chiral effective field theory formulated with 
a finite-range regulator provides a reliable and model-independent extrapolation to the physical world. 
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1 Introduction 

Quark nuclear physics describes our attempts to under- 
stand the structure of hadronic systems, including nuclei 
and dense matter, in terms of quarks and gluons - the 
fundamental degrees of freedom in QCD. It is impossi- 
ble in just a few pages to provide even a vague outline 
of the many exciting physics issues currently being ad- 
dressed in this field - from the possible phase transition 
to one or more quark-gluon phases at high temperature or 
density ^ , to suggestions of changes of hadron properties 
in-medium j3ini- Instead we shall concentrate on just one 
development which offers considerable insight into hadron 
structure from QCD itself, an approach which has led to 
surprisingly accurate comparisons between lattice QCD 
data and experiment as well as remarkable insights into 
how one might improve hadron models. 

As the time for calculations within lattice QCD j^ 
with dynamical fermions (including q — q creation and 
annihilation in the vacuum) scales as 'm~^'^ [S], current 
calculations have been limited to light quark masses 6-10 
times larger than the physical ones. With the next gen- 
eration of supercomputers, around 10 Teraflops, it should 
be possible to get as low as 2-3 times the physical quark 
mass, but to actually reach that goal on an acceptable vol- 
ume will require at least 500 Teraflops. This is probably 
10-20 years away. 

Since a major motivation for lattice QCD must be to 
unambiguously compare the calculations of hadron prop- 
erties with experiment, this is somewhat disappointing. 
The only remedy for the next decade at least is to find a 
way to extrapolate masses, form-factors, and so on, cal- 
culated at a range of masses considerably larger than the 
physical ones, to the chiral limit. In an eifort to avoid the- 
oretical bias this has usually been done through low-order 



polynomial fits as a function of quark mass. Unfortunately, 
as we discuss in Sec. 2, this is incorrect and can yield quite 
misleading results because of the Goldstone nature of the 
pion. 

The essential problem in performing calculations at 
realistic quark masses (of order 5 MeV) is the approxi- 
mate chiral symmetry of QCD. Goldstone's theorem tells 
us that chiral symmetry is dynamically broken and that 
the non-perturbative vacuum is highly non-trivial |B] , with 
massless Goldstone bosons in the limit niq -^ 0. For fi- 
nite quark mass these bosons are the three charge states 
of the pion with a mass m'^ oc niq. Although this re- 
sult strictly holds only for m^ near zero (the Gell-Mann- 
Oakes-Renner relation), lattice simulations show it is a 
good approximation for rn^ up to 1 GeV^ or so, and we 
shall use m^ here as a measure of the deviation from the 
chiral limit. 

On these very general grounds, one is therefore com- 
pelled to incorporate the non-analyticity into any extrap- 
olation procedure. The classical approach to this prob- 
lem is chiral perturbation theory (xpt), an effective field 
theory built upon the symmetries of QCD [7]. There is 
considerable evidence that the scale naturally associated 
with chiral symmetry breaking in QCD, yl^sB, is of order 
47r/^, or about 1 GeV. xPt then leads to an expansion in 
powers of mT^/A^sB and p/A^sB, with p a typical momen- 
tum scale for the process under consideration. At C'(p^), 
the corresponding effective Lagrangian has only a small 
number of unknown coefficients which can be determined 
from experiment. On the other hand, at 0{p^) there are 
more than 100 unknown parameters |8j, far too many to 
determine phenomenologically. 

While this situation seems formidable, the resolution 
is already in hand. We must first realize that the lattice 
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data obtained so far represents a wealth of information on 
the properties of hadrons within QCD itself over a range 
of quark masses. Just as the study of QCD as a function 
of Nc has taught us a great deal, so the behaviour as a 
function of m^ can yield considerable insight into hadronic 
physics. 

The first thing that stands out, once one views the data 
as a whole, is just how smoothly every hadron property 
behaves in the region of large quark mass. In fact, baryon 
masses behave like a + bniq, magnetic moments like (c + 
draq)~^, charge radii squared like (e + f mq)~^ and so 
on. Thus, if one defined a light "constituent quark mass" 
as M = Mq + crUq (with c ^ 1), one would find baryon 
masses proportional to M (times the number of u and d 
quarks), magnetic moments proportional to AI^^ and so 
on -just as in the constituent quark picture. There is little 
or no evidence for the rapid, non-linearity associated with 
the branch cuts created by Goldstone boson loops. 

Over the past few years we have come to a deep un- 
derstanding of why QCD exhibits these features. It will 
be helpful to summarise those conclusions here: 

— In the region of quark masses mq > 60 MeV or so (tott 
greater than typically 400-500 MeV) hadron properties 
are smooth, slowly varying functions of something like 
a constituent quark mass, M '^ Mo+cmq (with c '^ 1). 

— Indeed, Mjv ^ 3 M, Mp^^ ^ 2M and magnetic mo- 
ments behave like 1/Af . 

— As ruq decreases below 60 MeV or so, chiral symmetry 
leads to rapid, non-analytic variation, with 5Mn ^ 



m, 



3/2 



,5^lH 



rfif^l'^ and 



8 <r^ >ch"^ InrTiq. 
— Chiral quark models like the cloudy bag '^','10*11' pro- 
vide a natural explanation of this transition. The scale 
is basically set by the inverse size of the of the com- 
posite source, above which chiral loops are strongly 
suppressed. Below this scale the pion Compton wave- 
length is larger than the source and one begins to see 
rapid, non- analytic chiral corrections. 

These are remarkable results that will have profound 
consequences for our further exploration of hadron struc- 
ture within QCD as well as the analysis of the vast amount 
of data now being taken concerning unstable resonances. 
In terms of immediate results for the structure of the nu- 
cleon, we note that the careful incorporation of the cor- 
rect chiral behaviour of QCD into the extrapolation of its 
properties calculated on the lattice has produced: 

i) The most accurate values of the proton and neutron 

magnetic moments from lattice QCD |12| . 
ii) The most accurate value of the sigma commutator 

from lattice QCD fTSJ. 
iii) The most accurate values for the charge radii of the 

baryon octet from lattice QCD Il4j . 
iv) The most accurate values for the magnetic moments 

of the hyperons from lattice QCD 15,16 . 
v) Good agreement between the extrapolated moments of 

the non-singlet distribution, u — d, calculated in lattice 

QCD and the experimentally measured moments |17[ 

[El 



vi) The most accurate estimates of the low moments of 
the spin-dependent parton distribution functions at 
the physical quark mass from lattice QCD 19 . 

vii) An understanding of the failures of the assumption of 
universality of quark electromagnetic properties J2U| 
and an improved lattice estimate of the strangeness 
magnetic moment of the proton G%j 21 . 

Furthermore, this approach, together with the observed 
constituent-quark-like behaviour seen in the lattice data 
for TTiq > 50 MeV (as noted earlier), has suggested a novel 
way of modelling hadron structure |22II23| . 

Apart from the original publications, these develop- 
ments have been fairly widely reported at various confer- 
ences - e.g. see Refs. |23E1]- Here we focus particularly 
on the question of the extrapolation of hadron masses in 
order to clarify an issue of great current interest, namely 
the degree of model independence of the results obtained 
after a controlled chiral extrapolation. 



2 Effective Field Theory 

Chiral perturbation theory is a low-energy effective field 
theory for QCD. Low-energy properties of QCD can be 
expanded about the limit of vanishing momenta and quark 
mass. In relation to the extrapolation of lattice data, xpt 
provides a functional form applicable in the limit m^ — > 0. 

Goldstone boson loops give rise to specific corrections 
to baryon properties — most importantly, they give rise to 
non-analytic behaviour as a function of quark mass. The 
low-order, non-analytic contributions arise from the pole 
in the Goldstone boson propagator and hence are model- 
independent ^26j . Analytic variation of hadron properties 
is not constrained via the symmetry and hence expan- 
sions contain free parameters which must be determined 
by comparison with data. 

Effective field theory then tells us that the general ex- 
pansion of the nucleon mass about the SU(2) chiral limit 
is: 



niN = ao + a2 ml 



a4m* 



+(j^^{m^,A) + (JA7v{m-^,y^) + 



(1) 



where ctbtt is the self-energy arising from a Bn loop (B — 
N or A). The expansion has been written explicitly in this 
form to highlight that the theory is equivalently defined 
for an arbitrary regulator — see Ref. [27| for a complete 
discussion. 

The traditional approach within the literature is to use 
dimensional regularisation to evaluate the self-energy inte- 
grals. Under such a scheme the NNtt contribution simply 
becomes (TN7r(wT,yl) -^ clna^t^^ and the analytic terms, 
a„ m^ undergo an infinite renormalisation. The A con- 
tribution behaves similarly, producing a logarithm and 
one obtains a series expansion of the nucleon mass about 

771^ = 



mN = Co + C2 ml + C4 m^ 

-l-CLNA ml + CNLNA m^ In m^ -|- 



(2) 
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where the ai have been replaced by the renormalised (and 
finite) parameters q. 

It is not clear, a priori, that any such truncated expan- 
sion will be capable of reliably fitting lattice data. The first 
empirical indication of serious problems in this approach 
came with the realization that lattice data could not re- 
cover the model-independent coefficient, clna- Truncating 
the power series at the m^ term and allowing clna to vary 
as a free fit parameter, together with cq and C2, produced 
a value clna ^ —0.76 GeV~^ |2H|- This should be com- 
pared with the physical value of —5.6 GeV~^ — a factor 
of 8 larger! This tells us immediately that either there 
are serious convergence problems with the third order ex- 
pansion or lattice QCD is in error. Clearly most readers 
would opt for the first possibility and so do we. 

Even by retaining all terms as described in Eq. |5J| 
it is not clear that reliable extrapolation can be guaran- 
teed by fitting lattice data over a range of (heavy) quark 
masses. One point of issue is that it is derived in the limit 
771^ <^ A{= m/\ — TTiAf), whereas the lowest lattice data 
with dynamical fermions that one can expect in the next 
decade is perhaps 200-250 MeV - c.f. A = 292 MeV. It 
should be clear to those familiar with lattice simulations 
that even at this lightest pion mass the branch cut will 
not be observed due to the restricted phase space on a fi- 
nite volume lattice. Consequently all lattice data will still 
lie above A. Mathematically the region ttItt « Z\ is domi- 
nated by a square root branch cut which starts at m^ = A. 
Using dimensional regularisation this takes the form J2H|: 



6<?i 



257r2/.2 

A 



{A^ - TO^) 2 ln(Z\ + TO^ - V^2~ 



-{2A^ -3m^)lnm^ 



(3) 



for TOtt < A^ while for TOtt > A the first logarithm be- 
comes an arctangent. No serious attempt has been made 
to extend the formal expansion in Eq. (0) to incorporate 
this cut in an analysis of lattice data and, given the num- 
ber of parameters to be determined if one works to order 
m^, it is not likely that it will be done in the next decade. 

Even ignoring the Att cut for a short time, studies of 
the formal expansion of the N — >N tt ^ N self-energy inte- 
gral, ctntt, suggest that it has abysmal convergence proper- 
ties. Using a sharp, ultra-violet cut-off, Wright showed [3D] 
that the series diverged for m^ > 0.4 GeV. If one instead 
uses a dipole cut-off, which in view of the phenomenolog- 
ical shape of the nucleon's axial form-factor is much more 
realistic, it is worse - with the radius of convergence being 
around 0.25 GeV. We return to this in Sect. 3. 

The main issue of the convergence of this truncated 
series, Eq. ((2J), has its origin fixed in the formalism that 
it is derived from the general form of Eq. (QJ . The dimen- 
sionally regulated approach requires that the pion mass 
remain much lighter than every other mass scale involved 
in the problem. This requires that mT^/A^s's, ^ 1 and 
m^/Z\ <C 1. A further scale, as addressed in the intro- 
duction, is set by the physical extent of the source of the 
pion field. This scale, A ^ ^source' corresponds to the 



transition between rapid, non-linear variation and smooth, 
constituent-like quark mass behaviour. An alternative pro- 
cedure would be to regulate Eq. |^ with a finite A which 
physically corresponds to the source of the meson cloud 
having an extended structure. 

In summary, the low-energy effective field theory can 
be very useful in describing the quark mass behaviour of 
hadron properties. These powers have unfortunately been 
lost in the literature, where only a single type of regulator 
(i.e. dimensional) has been studied in detail. 



3 Accurate, Model-Independent Method of 
Chiral Extrapolation 

We now turn to the direct application of Eq. Q , written 
in a regulator independent form, to the extrapolation of 
lattice data. By studying a range of different regulators, 
both finite ranged and the dimensional approach, we can 
study the sensitivity to the form chosen. 

The additional analytic term, 014 m^, differs from pre- 
vious studies using a finite ranged regulator |5S]. Since 
one is working to non-analytic order m^ log m^ one is cer- 
tainly permitted freedom in a^ 27,. In practice, whether 
or not this can be reliably determined, together with A, 
will depend on the data available. 

The self- energies, ctbtt, are evaluated with a variety 
of regulators. Included here are both that of the trun- 
cated power series resulting from dimensional regularisa- 
tion and also those of finite range, including sharp cut-off, 
monopole, dipole and Gaussian. The coefficients of non- 
analytic terms are constrained to their phcnomcnological 
values, while the parameters a^ and A are then determined 
by fitting lattice data. 

It is again worth noting that, independent of regulator, 
precisely the same expansion of Eq. (21 will be obtained 
in the limit m^ <C A. Although the parameters ai will 
depend upon the choice of regulator, one can expand the 
self-energy terms to order m^ (or higher) to obtain the 
chiral coefficients at the appropriate order and compare 
with the truncated expansion of Eq. (O — see Ref. J2Z1 
for a full discussion of this issue. 

The results of our fits to lattice data studied with 
a range of regulators are shown in Fig. ^ The short- 
dash curve shows the fit obtained using the dimension- 
ally regulated form of Eq. (21 extended to a further term 
in the expansion, aQm%. This additional term is neces- 
sary for a more reasonable fit due to the large NLNA 
contribution. The long dash curve is the result of using a 
sharp cut-off regulator, while the following three fits (solid 
curves), indistinguishable on this plot correspond to 
the monopole, dipole and Gaussian regulated forms. 

The extrapolation of lattice data is clearly independent 
of the choice of finite, ultra-violet regulator. Knowing this 
we can examine the range of convergence of a truncated 
power series. Selecting the dipole form and neglecting the 
Att contribution to the nucleon self energy one can obtain 
a closed analytic expression. An expansion in powers of the 
pion mass is shown in Fig. (21 Here we see that convergence 
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Fig. 1. Fits to lattice data j3i. for five difi'erent ultra-violet 
regulators. 




0.4 0.6 

Fig. 2. The fit to the lattice data using the dipole regulator. 
The dashed curves show power series expansions of this fit to 
successive orders in m„ for m^ -^ m%. 



of the expansion to order m^ breaks down above ttItt ^ 
250 MeV. 

For the reasons outlined, it is essential that the self- 
energies are evaluated using some ultra-violet regulator — 
a sharp cut-off or a dipole form, for example. Whatever 
is chosen does not effect the non-analytic structure which 
is guaranteed correct. The branch points at ttItt equals 
zero and A are incorporated naturally. The use of a finite 
regulator then automatically produces the transition scale 
associated with the physical extent of the meson source. 

The essential point is that studies of the nucleon (c.f. 
Ref. ISni and Fig.[TJ, the A (c.f. Fig. 4 of Leinweber et al. 
|28j ') and the p meson '32' suggest that this procedure will 
result in little or no model dependence in the extrapolation 
to the physical pion mass once there is accurate lattice data 
for niTr ~ 0.3 GeV or less. Physically this is possible be- 
cause the self-energy loops are rapidly suppressed in the 
region ttItt > 0.4 GeV. Thus, an extrapolation based on 
Eq. (^ formulated with the selection of a long-distance 
regulator allows one to respect all the chiral constraints, 
keep the number of fitting parameters low and yield es- 



sentially model-independent results at the physical pion 
mass. No other approach can do this. 



4 Possible Connection to QQCD 

Although quenched QCD (QQCD) is an unphysical the- 
ory, it provides an alternative avenue for enhancing our 
understanding of chiral extrapolation. Multi-mass tech- 
niques allow a dense set of quark masses to be simulated 
with relative computational ease. Together with recent ad- 
vances in numerical techniques, which allow simulations 
to be performed at light quarks masses ^33! , one will be 
able to very accurately determine the quark mass depen- 
dence of quenched simulations within the light quark mass 
regime. 

The study of baryon spectroscopy in quenched lat- 
tice QCD has recently made great progress. We have al- 
ready noted that the lattice data behaves like a constituent 
quark model for quark masses above 50-60 MeV because 
Goldstone boson loops are strongly suppressed in this re- 
gion. This not only provides a very natural explanation 
of the similarity of quenched and full data in this region 
but it also suggests a much more ambitious approach to 
hadron spectra. It suggests that one might remove the 
small effects of Goldstone boson loops in QQCD (includ- 
ing the rj') and then estimate the hadron masses in full 
QCD by introducing the Goldstone loops which yield the 
LNA and NLNA behaviour in full QCD. 

As a first test of this idea, Young et al. j33] recently 
analysed the MILC data [31] for the N and A, using Eq.© 
(with a4 = 0) for full QCD and the appropriate general- 
ization for QQCD - i.e. using quenched pion couplings as 
well as the single- and double- "hairpin" 77' loops [SEIEZI- 
The results illustrated in Fig. |3| are remarkable. The val- 
ues of ao and 0:2 for the N (or the A) obtained in QQCD 
agree within statistical errors with those obtained in full 
QCD. Certainly this result (unlike the result for the ex- 
trapolation of individual hadron masses as noted above) 
is somewhat dependent on the shape of the ultra-violet 
cut-off chosen - although the extent of that is yet to be 
studied in detail. Nevertheless, given that the study in- 
volved the phenomenologically favoured dipole form, it is 
a remarkable result and merits further investigation. 



5 Conclusion 

At the present time we have a wonderful conjunction of 
opportunities. Modern accelerator facilities are providing 
data of unprecedented precision over a tremendous kine- 
matic range at the same time as numerical simulations of 
lattice QCD are delivering results of impressive accuracy. 
It is therefore timely to ask how to use these advances to 
develop a new and deeper understanding of hadron struc- 
ture and dynamics. 

We have demonstrated that the use of chiral effective 
field theory can provide accurate extrapolation formulae. 
In particular, we have shown that the extrapolation of the 
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Fig. 3. Fits to both quenched (open triangles) and dynamical 
(closed triangles) lattice data |35| using a dipole regulator 34 . 



nucleon mass exhibits minimal model-dependence in the 
choice of finite-ranged, ultra-violet regulator. 

In combination with the very successful techniques for 
chiral extrapolation, lattice QCD will finally yield accu- 
rate data on the consequences of non-perturbative QCD. 
Furthermore, the physical insights obtained from the study 
of hadron properties as a function of quark mass will guide 
the development of new quark models and hence a much 
more realistic picture of hadron structure. 
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